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Physical configurations corresponding to topological excitations present in the XY 
limit of a quantum spin | Heisenberg ferromagnet, are investigated on a two dimensional 
square lattice . Quantum vortices (anti-vortices) are constructed in terms of coherent spin 
field components and the crucial role of the associated Wess-Zumino term is highlighted . 
It is shown that this term can identify a large class of vortices (anti- vortices). In particular 
the excitations with odd topological charge belonging to this class , are found to exhibit a 
self-similar pattern regarding the internal charge distribution . Our formalism is distinctly 
different from the conventional approach for the construction of quantum vortices ( anti- 
vortices ). 
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Introduction 

The physical existence of topological excitations in quantum ferromagnetic spin sys- 
tems in two dimension can be investigated by making use of coherent state formulation 
1,2,3 rp^g XY-limit of Heisenberg ferromagnet also belongs to this class of systems . It has 
been shown that in the case of a Heisenberg antiferromagnetic chain, the Wess-Zumino ( 
WZ) term in the effective action is genuinely a topological term l . Moreover in the long 
wavelength limit the explicit expression of this term is similar to that of the topological 
term for a nonlinear sigma model. Extending this approach further, it was shown by us 
3 that it is possible to express the WZ term as a topological term, also in the case of a 
ferromagnetic chain in the long wavelength limit. We generalized this to the case of a 2D 
ferromagnet 3 ' 4 to demonstrate that in the above limit we could indeed get an expression 
from the WZ term , indicating the possibility of topological excitations. 

Thus at this stage of analysis it is necessary to examine the physical configurations of 
these excitations in terms of coherent fields as well as the consequence of the WZ term, 
explicitly on a two dimensional square lattice. At the moment, we do this for the XY-limit 
of the 2D Heisenberg spin | ferromagnet. We show that the WZ term really identifies 
a large class of the topological excitations . Furthermore we demonstrate that this term 
can clearly differentiate between vortices ( and antivortices ) with different charges within 
this class. It is mentionworthy that our entire analysis is valid for all temperatures . 
Mathematical Formulation 

The quantum Euclidean action S^fn] for the spin coherent fields n(t) can be written as 2,3 

S E [n] = -isS W z[n] + — / dt[d t n(t)] 2 + / dtH(n) (1) 

4 jo Jo 

where s is the magnitude of the spin (s = = in the present case ) and 

ff(n) = <n|fT(S)|n) (2) 
H(S) being the spin Hamiltonian in representation V. In the limit 5t — > the term 

2 



^ Jq dt[d t n(t)] 2 on the right hand side of Eqn.(l) vanishes and henceforth we omit this 
term in the expression for the action. The Wess-Zumino term Swz is given by 2 

rfi rl 

S wz [n] = dt drn.it, r) ■ d t n(t, r) A d T n(t, r) (3) 
Jo Jo 

with n(t, 0) = n(t), n(t, 1) = n , and n(0, r) = n(/3, r), t G [0, 0], r G [0, 1] , /? being the 
usual inverse temperature. 

The expression in equation(3) is the area of the cap bounded by the trajectory T parametrized 
by n(t) [=(ni(t), n 2 (t), n 3 (t))] on the sphere: 

n ■ n = 1 (4) 

Here |n) is the spin coherent state on a single lattice point as defined in Refs. 1-3 and 
all the expressions above are on a single lattice point. We calculate the difference {5 Swz) 
between the WZ terms on a pair of neighbouring lattice points 2 and express the WZ term 
for the whole lattice as 

2N 

s wz = J2 S wzHiaJa)} 

2N 1 l 

= H b^Kt' ~ l ) a i3 a )] + 2 S wz[n{ia, (j ~ l)a)] 

1 1 

+ -5 x S W z[n-(ia, ja)} + -5 y Swz[n-(ia, ja)]} (5) 



where 



S x S W z[n(T)] = S W z[n(iaJa)]-S wz [n{(i-l)a,ja)] 
dt[8 x n-(nAd t n)](r) 



o 



5 y Swz[n(r)] = S wz [n(ia, ja)] - S wz [n{ia, U - 1)«)] 

= [ dt[5 y n-{nAd t n)](r) (6) 

J 

r = (ia,ja) , a being the lattice spacing and i,j = 1, 2, , 2N. Now the spin Hamiltonian 

corresponding to anisotropic Heisenberg ferromagnet of XX Z type is given by 

H(S) = -<?£S(r) ■ S(p/) - g\ z E S m (t)S m (ti) (7) 

(r,r) ( r i r ') 



with g > and < X z < 1 , r,r/ run over the lattice, and (r, r/) signifies nearest neighbours 
and S = (S, S z ). 

It follows from Equations (2) and (7) that the spin Hamiltonian in terms of coherent fields 
is given by 

H(n) = -g n(ia,ja)-h(ra,sa)-X z n z (ia,ja)n z (ra,sa) (8) 

((Vi).(r,«)) ((i,j),(r,s)) 

Now the quantum action (Euclidean)for the two dimensional anisotropic ferromagnet is 
given by 

rP 

S E [n] = -isYSwz[n(ia,ja)]+ dt[-gs 2 J2 n(z, j) ■ n(r, s) 

hi ° (W),(r,«)> 

-g\ z s 2 Y n z (ia,ja)n z (ra,sa)] (9) 

<(i,i).( r > s )> 

where s — \ since we are considering extreme quantum case. 

In order to evaluate WZ-term [ i.e.,S^ z ] on a vortex plaquette, we need to know the 

expression of the WZ-term on the lattice with the help of Eqs.(5) & (6). 

Now 

S x n(ia,ja) = n(ia, ja) - n((i - l)a, ja) 

5 y ri(ia,ja) = n(ia,ja) — n(ia, (j — l)a) (10) 
From Eqns. (5) , (6) and (10) we get 

2S wz [n(ia,ja)] = S wz [n((i - l)a, ja)] + S wz [n(ia, (j - l)a)] 

+ / dt[n(ia, ja) — n((i — l)a, ja)} ■ (n A d t n)(ia, ja) 
Jo 

rH 

+ / dt[n(ia, ja) — n(ia, (j — l)a)] ■ (n A d t n)(ia, ja) 
Jo 

= S wz [n((i - l)a, ja)] + S wz [n(ia, (j - I) a)] 
rP 

— I dtn{{i — l)a,ja) ■ (n A d t n)(ia, ja) 
Jo 

rP 

— / dtn(ia, (j — l)a) ■ (n A d t n)(ia, ja) (11) 
Jo 

The time derivative of the coherent spin field, <9jn in the above equation can be expressed 
in terms of coherent fields and their spatial derivatives through equations of motion. This 



we obtain from the action given by Eq.(9) in the continuum limit as an approximation . 
Thus in this limit we have, 

S tot = S E [n] + J d 2 x^A(x,t)(n 2 (x,t) - 1) (12) 

where the first action in the right hand side of the above equation is the continuum 
version of the action given by Eq.(9) . A(x, t) in the above equation is an auxiliary field 
playing the role of a multiplier . The equations which follow from the above action are 
the following 

d t rn = (nAV 2 n) 1 -(l-A,)n 2 V 2 n 3 



2ga? 
i 



d t n 2 = (nAV 2 n) 2 + (l-A>iV 2 n 3 



2ga? 

% -d t n 3 = (nAV 2 n) 3 (13) 



2ga 

Now we substitute the expressions for c^n from the above equations, into Eq. (11) with 
discretized version of the derivatives on the right hand side of Eq. (13) . This enables 
us to obtain an expression for Sffiz corresponding to the XXZ Heisenberg ferromagnet. 
The expressions of the derivatives to be used on the two dimensional lattice are as follows 
[assuming the common factor ^2 to be equal to 1 for simplicity] :- 

d x ii(iaja) = -[n(ia,ja)-n((i-l)a,ja)] 
a 

d y n(ia,ja) = -[n(ia, ja) - n(io, (j - l)o)] (14) 
V 2 n(ia,ja) = — [2n(za, ja) — 2n((z — l)a, ja) — 2n(za, (j — l)a) 



a 



+n((i - 2)o, jo) + n(ia, (j - 2)a)] (15) 

To keep the calculations simple but reasonable , we retain the intra-plaquette contributions 
by imposing a "local periodic boundary condition" (local PBC) viz., 

n(ia,ja) = n(ia + (Q + l)ia,ja) 

n(ia,ja) = n(ia, (Q + l)ja + ja) (16) 



where Q is the "topological charge" of the plaquette under consideration. We use the 
above condition in the following section where we consider vortices with topological 
charges 1,2,3 etc. Henceforth we denote a vortex with topological charge Q ( a con- 
figuration on the lattice in which the spin undergoes a 2tc ■ Q rotation as we go around 
the plaquette once ) as Q-vortex . 



Calculations and Results 
: Analysis of 1- vortex : 

Eq.(ll) shows that the last two terms are the contributions from S-wziia-, 3 a ) to the 




Fig. la Fig. lb 

Figure 1: (a) 1-vortex , (b) 1-anti- vortex 

vertices 2 , 3 and 4 with coordinates (ia, (j — l)a) , (ia,ja) and ((i — l)a,ja) respectively 
, belonging to the 1-vortex [see Fig.l]. Other contributions to the plaquette can be 
obtained by considering Swz at different lattice points. In other words , for a vortex ( 
anti-vortex ) plaquette having vertices 1,2,3 and 4 with coordinates ((z — l)a, (j — l)a), 
(ia, (j — l)a), (ia,ja) and ((i — l)a,ja) respectively, Swz associates the kind of terms, 
similar to the ones discussed above , to a 1-vortex ( 1-antivortex ) [seeFig.l]. It is easy to 
check that the contribution from Swz to a single plaquette (with vertices ((«— l)a, (j—l)a), 
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(ia, (j — l)a), (ia,ja) and ((z — l)a, ja) are given by (see Fig.l ) 

V\ = C{ dt[n((i — l)a,ja) ■ (n A d t n)(ia, ja) 
Jo 

+n(ia, (j — l)a) ■ (n A d t n)(ia, ja) 

+n((i - l)a, (j - l)o) ■ (n A d t n)(ia, (j - l)a) 

+n((i - l)a, (j - l)o) • (n A <9 t n)((z - l)a, ja)]} (17) 

where the proportionality constant C depends on the size of the lattice . Henceforth we 
drop C from all the expressions for convenience. Since we will be interested in the extreme 
XY- anisotropic limit we we assume n 3 = sin e on each lattice point where e is very small. 
In this equation one of the components n± , n 2 is given by 1 — 5 so that the other is given 
by ±V cos 2 e — (1 — 5) 2 , 5 being small . In Figure-1 we have a quantum vortex in which 
the horizontal arrow — > at a vertex implies n\ with value 1 — 5 and the vertical arrow j 
implies n 2 with value 1 — 5. Also in this figure the horizontal arrow <— at a vertex implies 
ni with value —(1 — 5) and the vertical arrow j implies n 2 with value —(1 — 5). 
In this connection let us point out that usually in a two dimensional vortex corresponding 
to spin I quantum spin model the arrows — > and f signify eigenstates of S x and S"y with 
eigenvalues +| respectively 5 . <— and | are eigenstates of these operators with eigenvalues 
— I respectively. In our quantum formulation we have the spin coherent field components 
n\ and n 2 ( satisfying the constraint given by Eqn.(4) at each lattice point ) forming the 
vortex (anti-vortex) with 77.3 = sine at each lattice point. Our picture is that of a flattened 
"meron configuration" 6,7 mimicking a conventional vortex or anti- vortex . Therefore n\ or 
n 2 cannot be exactly equal to ±1 . We choose them to be ±(1 — 5), where 5 is a function of 
e with 5 — > as e — ► 0. Moreover if we want to identify the spin state | — >•) (eigenstate 
of the operator S x with eigenvalue +|) with the coherent state cos 1 1|) + (e~^)sin || — |) 
at a vertex [Fig.l] we have n± = l,n 2 = and 77,3 = at that vertex. Similarly for | <— ) 
we have ni = — 1, 72,2 = and n 3 = ; for | |) we have ni = 0, n 2 = 1 arifi 77.3 = ; for 
I J,) rii = 0, n 2 — — 1 and n 3 = , in case we want to identify the respective spin states 
with the coherent spin states. This does not agree with nonzero magnitude of 77.3 , given 
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by n 3 = sine . Thus we choose n\ or n 2 to be ±(1 — 5) to preserve the constraint given 
by Eqn.(A). 

Now it can be easily verified from Eqn.(15) that V 2 n 3 = at all lattice points. Then using 
Eqns. (11) and (13) we can easily show that the contribution of Sffiz t° the plaquette in 
Fig.l is 



+ [-A(n 2 2 (3) + sin 2 e) + 2sin 2 e ri!(3)]VV(3) 

+ [ni(l)n 2 (2)ni(2) - n 2 (l)[sin 2 e + m 2 (2)] + sin 2 e n 2 (2)]V 2 ri 2 (2) 

+ [-rii(l)(sin 2 e + n 2 2 (2)) + n 2 (l)n 2 (2)m(2) + sin 2 e m(2)]V 2 ni(2) 

+ [m(l)n 2 ni(4) - n 2 (l)(sin 2 e + rii 2 (4)) + sin 2 e n 2 (4)] V 2 n 2 (4) 

+ [-r2 1 (l)(sin 2 e + n 2 2 (4)) +n 2 (l)n 1 n 2 (4) + sin 2 e m(4)]VV(4)} (18) 



Using Eqns. (15) and (16) the derivatives on the right hand side of .Eg. (18) are given by: 



V 2 n 2 (3) = 4n 2 (3) 

V\(3) = 4m(3)-2« 

V 2 n 2 (2) = 4n 2 (2)-2/x 

V 2 n 2 (2) = 4th (2) 

V 2 n 2 (4) = 4n 2 (4)-2 / u 

VV(4) = 4m(4) (19) 



Notice that while obtaining Eq. (18) we have assumed the local PBC [Eqn. (16)} so that 
the derivative V 2 n at a point on the plaquette [Fig.l] is written in terms of the fields 
n at points on the same plaquette. In Eq. (18) we have also made use of the fact that 




where /x and n are given by 



n 2 (l) + n 2 (3) 



ni(2) + m(4) 



(20) 



S 



7ii(l)+rii(3) = (1 -5) - (1-6) = and n 2 (2)+n 2 (3) = (1 - 5) - (1 - 5) = in [Fig!]. 
Now we subtitute Eq.(l&) into Eq. (17) and write SVzi as A + S. It is interesting to note 
that A remains the same if we go from vortex to antivortex by changing n 2 (2) and n 2 (4) 
in Fig A to — n 2 (2) and — n 2 (4) respectively whereas goes over to —B . Thus Sy/Zx 
takes the form A — B for antivortex. We have the following expression for Swz 1 for a 
1 -vortex . 

Swz x = A + B 

A = / dt{8sin 2 e cos 2 e 
Jo 

+2/in 2 (l)[2sin 2 e + m 2 (2) + m 2 (4)] + 4sin 2 e[n 2 2 (2) + n 2 2 (4)]} 

B = -4 / dt /J, m(l)ni(2)n 2 (3) -n 2 (2) (21) 
Jo 

where we have used the symmetry property that n 2 (2) = — n 2 (4) and further the constraint 
n 2 = 1 at every site. It should be noted that the entire scheme is meaningful only if B 
is nonzero. The case B = corresponds to the configurations in Fig. 1 with n 2 (l) = 
ni(2) = n 2 (3) = ni(4) = . We consider only those configurations in Fig. 1 in which 
none of the above equations is satisfied . The WZ term distinguishes between vortex or 
antivortex for such configurations only , corresponding to the above case . We will see 
more such B = configurations corresponding to higher topological charges , for which 
WZ term is not sensitive. 

: Analysis of 2- vortex : 
For a typical 2- vortex we refer to Fig. 2 . We first calculate the contribution of the WZ- 
term on such a plaquette by simply adding the contributions of WZ from each of the 
individual elementary plaquettes (subvortices) viz., a , b , c and d, each carrying topolog- 
ical charge +1 , belonging to the above 2- vortex. This we denote by [a + b + c + d] FREE . 
We also calculate the contributions of the WZ-term from the above four subvortices as 
they are glued in a natural way to form the above mentioned 2-vortex [ see in Fig.2] . We 
implement this by taking the contribution of the WZ-term with a factor of \ by taking 
into consideration the common bonds occuring between the pairs of adjacent subvortices. 

9 
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Figure 2: 2- vortex 

We denote this total contribution by [a + b + c + d] Gi!mD . Let us now write down the 
formal expressions of [a + b + c + d] FflEE and [a + b + c + d\ GLUED . 

r/3 

[a + 6 + c + c/] FR£;£; = / dt{n(la) • (n A <%n) (2a) + n(2a) • (n A d t n) (3a) 

«/ 

+n(4a) ■ (n A 5 t n)(3a) + n(la) • (n A d t n)(Aa) 



+n(16) 
+n(46) • 
+n(lc) • 
+n(4c) • 
+n(ld) • 
+n(4d) • 



n(16) • (n A d t n)(2b) + n(26) • (n A d t n)(3b) 
n A d t n)(3b) + n(16) • (n A <9 t n)(46) 
n A d t n)(2c) + n(2c) • (n A d t n)(3c) 
n A 9 t n)(3c) + n(lc) • (n A 9*n)(4c) 
n A d t n)(2d) + n(2d) ■ (n A d t n)(3d) 
n A d t n)(3d) + n(ld) • (n A ftn)(4d)} 

r/3 1 
[a + b + c + d] GLUED = / dt{n(la) • (n A d t n)(2a) + -n(2a) • (n A dtn) (3a) 

JO Zi 

+ ^n(4a) • (n A d t n) (3a) + n(la) ■ (n A d t n) (4a) 
+n(16) • (n A d t u)(2b) + n(26) ■ (n A d t n)(3b) 



(22) 
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+ in(46) • (n A d t n)(3b) + -n(16) • (n A d t n)(4b) 
+ 2 n ( lc ) • ( n A 9 t n)(2c) + n(2c) • (n A d t n)(3c) 
+n(4c) • (nAd,n)(3c) + ^n(lc) ■ (n A 0*11) (4c) 
+-n(ld) ■ (n A d t n)(2d) + -n(2d) ■ (n A $n)(3d) 
+n(4d) • (nAa t n)(3d) + n(ld) • (n A flfcn) (4d) } (23) 

With the help of Eqs. (13) and (15), we now determine the following quantity R: 

R = RELEVANT PART{[a + b + c + d] FREE - 2[a + b + c + d] GLUED } (24) 

The RELEVANT PART contains the terms which are linear in n 2 (2), 77-2(4), n 2 (6) and 
712(8). They are called RELEVANT because if we flip the directions of 712(2), Tt 2 (4), 7i 2 (6) 
and 72,2(8) , we go from the state of vortex to its corresponding anti-vortex state [see Fig. 2] 
.The quantity R has been chosen in a way so as to represent the degree of gluing amongst 
the constituent elementary plaquettes , while the composite object undergoes topological 
charge reversal. Thus the vanishing of R provides a minimal condition for the lattice spin 
configuration on a plaquette to attain a topological character . 
Now the quantity R defined in Eq. (24) is given by 

R = dt{n 2 (2)[ ni (4) -ni(6)][-2ni(2)n 2 (l) + n 1 (2)n 1 (3) 
Jo 

+4m(l)n 2 (3) + 2n 2 (0)n 1 (l) - 7i 2 (5)n 1 (l)]} (25) 

where we have assumed for simplicity ni(2) = ni(8). Thus the right hand side of the 
above equation vanishes if the following two conditions are satisfied 

m(2) = 711(8) 

m(4) = m(6) (26) 

These conditions can easily be satisfied for a large class of spin configurations possessing 

topological charge 2. 

Henceforth we consider the equation 

R = (27) 
11 



as the criterion for a spin configuration corresponding to a Q-vortex, to be a topological 
excitation , satisfying the internal consistency . 
For Q=3 [see Fig. 3], for example ,we have 

R = RELEVANT PART {[a + b+ .... + h] FREE - 2 [a + b + .... + h} GLUED } (28) 

with Eqn. (27) being the criterion for the 3-vortex to be a topological excitation . 

Let us mention that in the present case for Q=2, there are configurations correspond- 
ing to the vanishing of the relevant term in the glued vortex i.e, B — , which fulfil 
the condition R = as well in Eqn. (27). Our scheme is however applicable to only 
R = configurations having finite B . The fact that the WZ term is sensitive only to a 
particular class of configurations , has its origin in the presence of a non-topological part 
along with the topological one [see Ref. 3] , in it. 

:Analysis of 3-vortex : 

We directly refer to the Fig. 3. We obtain the cofugurations corresponding to the 
following minimal conditions satisfying Eqn. (27) 

m(2) = m(i2) 

m(6) = m(8) 

m(4) = m(io) 

m(B) = 

n 2 (l) = n 2 (3) 

n 2 (5) = n 2 (9) 

n 2 (3) = n 2 (ll) 

n 2 (C) = n 2 (D) 

Here also we do not consider the configurations corresponding to B 
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(29) 

in the glued 
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Figure 3: 3- vortex 
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vortex. Studying the configuration in Fig. 3 we discover the following identity for a 
topological Q- vortex ,with odd Q ( Q ^ 1 ), which describes the topological charge 
ditribution inside the vortex consistently. 

1 



Q = -[Q 2 -(Q-2) 2 ] + l 



(30) 



The quantity (Q — 2) is the absolute value of the effective charge ( negative ) of the core 
which is an anti-vortex in this case. Notice that the above equation holds good even within 
the core, i.e., when Q is replaced by (Q — 2), and thereby depicts a self similar pattern. 
:Analysis of 4-vortex : 
In the case of 4-vortex Fig. 4 we have : 

R = RELEVANT PART {[a + b+ + k + l] FREE - 2[a + b + + k + l] GLUED } 

(31) 

and minimal solutions of Eqn. (27) as in the following : 



ni(even points) = ri2(odd points) = 
13 



(32) 
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Figure 4: 4-vortex 

where even and odd points imply the vertices on the boundary of the 4-vortex with even 
and odd integer labels. This solution [Eqn. (32)] also corresponds to vanishing of the term 
B = for the glued vortex. Moreover, the self similarity pattern, as discussed in Eqn. (30) 
, breaks down in this case. 

The above two features continue to hold for vortices ( anti- vortices ) of higher even- valued 
topological charges . 

Conclusion and Discussion 
i) Our calculations and analysis with 2D spin | ferromagnetic quantum XY model clearly 
bring out distingushing features between the even and odd charge sectors. The charge 
distribution equation [Eqn. (30)] involving the full spin configuration, is obeyed by the 
odd charge excitation in totality. The even charge sector however is not governed by this 
equation. 

ii) Our work has established the role of WZ term as a topological charge measuring 
quantity , obtained from microscopic theory. Furthermore, this term is able to test var- 
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ious internal consistencies and conservation conditions involving the topological charge 
distribution for a large class of composite excitations . Thus our approach is more pow- 
erful than that based on the heuristic operators suggested for determining the charges of 
topological excitations 5 ' 8 . 

iii) The restriction on the class of excitations identifiable by WZ term on a lattice , 
arises from the presence of a non-topological contribution , besides the usual topological 
one 3 . 

iv) Our future plan includes the generalization of our approach to the case of finite 
A 2 to achieve a physical realization of meronic type of excitations in the spin models . 
Furthermore making use of these results , we aim to evaluate the static and dynamic 
spin correlations for two dimensional spin | anisotropic quantum Heisenberg ferromagnet 
at any temperature . This can also throw some light on the possible phenomenological 
scenario of quantum Kosterlitz - Thouless transition 3 ' 9 . 

To conclude, this study of topological spin excitations on the 2D-lattice will undoubtedly 
play an important role in the understanding of thermodynamics of low dimensional fer- 
romagnets. 

Acknowledgement : One of the authors (SKP) would like to thank M.G.Mustafa, Ra- 
jarshi Roy and Purnendu Chakraborti for their valuable help in the preparation of the 
figures. 

We dedicate this piece of work to the memory of Late Professor Chanchal Kumar 
Majumdar. 
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